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Abstract 

We report preliminary results of a Dalitz plot analysis of three-body charmless B° — > K^n+ir^ 
decays. The analysis is performed with a data sample that contains 386 million BB pairs col- 
lected near the T(4S*) resonance with the Belle detector at the KEKB asymmetric energy e + e" 
collider. Measurements of branching fractions for quasi-two-body decays B° — > p(770)°K°, 
B° -» f (980)K°, B° -> K*(892)+ir-, B° -> iT"(1430) + 7r-, and upper limits on several other 
quasi-two-body decay modes are reported. 

PACS numbers: 13.20.He, 13.25.Hw, 13.30.Eg, 14.40.Nd 
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INTRODUCTION 



First results of amplitude analyses of B meson decays to a several three-body charmless 
hadronic final states have been reported recently: B + — > K + K + K~ [1], B + — > K + ii + 7i~ 0, 
HI H], B° — > K + 7i^7i° 0] and B + — ► 7r + 7r + 7r~ j^. Branching fractions for a number of 
quasi-two-body decays have been measured with some of them being observed for the first 
time. 

In this paper we present preliminary results of a Dalitz plot analysis of neutral B meson 
decay to the K^-n +r n^ three-body charmless final state. The analysis is based on a 357 fb _1 
data sample containing 386 million BB pairs, collected with the Belle detector operating at 
the KEKB asymmetric-energy e + e~ collider with a center-of-mass (cm.) energy at the 
T(4S) resonance (on- resonance data). The beam energies are 3.5 GeV for positrons and 8.0 
GeV for electrons. For the study of the e + e~ — > qq continuum background, we use data 
taken 60 MeV below the T(4S") resonance (off-resonance data). 

THE BELLE DETECTOR 

The Belle detector Q is a large-solid-angle magnetic spectrometer based on a 1.5 T 
superconducting solenoid magnet. Charged particle tracking is provided by a silicon vertex 
detector and a 50-layer central drift chamber (CDC) that surround the interaction point. 
Two inner detector configurations were used. A 2.0 cm beampipe and a 3-layer silicon vertex 
detector was used for the first sample of 152 million BB pairs, while a 1.5 cm beampipe, a 
4-layer silicon detector and a small-cell inner drift chamber were used to record the remaining 
234 million BB pairs 0]. The charged particle acceptance covers laboratory polar angles 
between 9 = 17° and 150°, corresponding to about 92% of the total solid angle in the cm. 
frame. The momentum resolution is determined from cosmic rays and e + e~ — > fi + n~ events 
to be cr pt /p t = (0.30 © 0.19pi)%, where p t is the transverse momentum in GeV/c. 

Charged hadron identification is provided by dE/dx measurements in the CDC, an ar- 
ray of 1188 aerogel Cerenkov counters (ACC), and a barrel-like array of 128 time-of-flight 
scintillation counters (TOF); information from the three subdetectors is combined to form 
a single likelihood ratio, which is then used in kaon and pion selection. Electromagnetic 
showering particles are detected in an array of 8736 CsI(Tl) crystals (ECL) that covers the 
same solid angle as the charged particle tracking system. The energy resolution for elec- 
tromagnetic showers is a E /E = (1.3 © 0.07/E © 0.8/£ 1/4 )%, where E is in GeV. Electron 
identification in Belle is based on a combination of dE/dx measurements in the CDC, the 
response of the ACC, and the position, shape and total energy deposition (i.e., E/p) of the 
shower detected in the ECL. The electron identification efficiency is greater than 92% for 
tracks with p^ > 1.0 GeV/c and the hadron misidentification probability is below 0.3%. 
The magnetic field is returned via an iron yoke that is instrumented to detect muons and 
K\ mesons. We use a GEANT-based Monte Carlo (MC) simulation to model the response 
of the detector and determine its acceptance 

EVENT RECONSTRUCTION 

Charged tracks are selected with a set of track quality requirements based on the number 
of CDC hits and on the distances of closest approach to the interaction point. We also 
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require that the track momenta transverse to the beam be greater than 0.1 GeV/c to reduce 
the low momentum combinatorial background. Charged tracks that are positively identified 
as kaons, protons or electrons are excluded. Since the muon identification efficiency and fake 
rate vary significantly with the track momentum, we do not veto muons to avoid additional 
systematic errors. 

We identify B candidates using two variables: the difference AE between the total re- 
constructed energy of a three-body combination and the nominal beam energy in the cm. 
frame and the beam constrained mass M hc . AE is calculated as AE = E B — -^bcam = 



\J2i y c2 P 2 + c 4 mfj — -Ebcam! where the summation is over all particles from a B candidate; 

and Pi and rrii are their cm. three- momenta and masses, respectively. Since there are no 
7r°'s or photons in the final state, the AE width (with typical value of 15 MeV) is governed 
by the track momentum resolution. The beam energy spread is about 3 MeV and gives a 
negligible contribution to the total AE width. The signal AE shape is parametrized by a 
sum of two Gaussian functions with a common mean. The AE shape for the qq background 
is parametrized by a linear function. The beam constrained mass variable Mt, c is equivalent 
to the B invariant mass with the measured B candidate energy Eb replaced by the beam 



energy E* cam . M hc = ±^Eg mn -c*P B = ± y/Eg^ - c 2 (E 4 p) 2 , where P B is the B can- 



didate momentum in the cm. frame. The average B meson momentum in the cm. frame is 
about 0.34 GeV/c which is much smaller than its total energy. Thus, the uncertainty in the 
measured P# gives a small contribution to the M bc width, which is dominated by the beam 
energy spread. The M bc width is about 3 MeV/c 2 and well described by a single Gaussian 
function. The M bc width in general, does not depend on the final state (unless photons are 
included in the reconstructed final state). 

BACKGROUND SUPPRESSION 

There are two sources of the background: the dominant one is due to e + e~ — > qq 
(q = u,d,s and c quarks) continuum events that have a cross-section about three times 
larger than that for the e + e~ — > T(AS) — > BB; the other one originates from other B 
meson decays. The background from continuum events is suppressed using variables that 
characterize the event topology. Since the two B mesons produced from an T(4S*) decay are 
nearly at rest in the cm. frame, their decay products are uncorrelated and the event tends 
to be spherical. In contrast, hadrons from continuum qq events tend to exhibit a two-jet 
structure. We use 6? t h r , which is the angle between the thrust axis of the B candidate and 
that of the rest of the event, to discriminate between the two cases. The distribution of 
(cos^thrl, is strongly peaked near Icos^thrl = 1-0 for qq events and is nearly flat for BB 
events. We require |cos# t hr| < 0.80 eliminating about 83% of the continuum background 
and retaining 79% of the signal events. For further suppression of the continuum background 
we use a Fisher discriminant formed from 11 variables: nine variables that characterize the 
angular distribution of the momentum flow in the event with respect to the B candidate 
thrust axis, the angle of the B candidate thrust axis with respect to the beam axis, and 
the angle between the B candidate momentum and the beam axis. The discriminant, is 
the linear combination of the input variables that maximizes the separation between signal 
and background. The coefficients are determined using off-resonance data and a large set 
of signal MC events. Use of the Fisher discriminant rejects about 89% of the remaining 
continuum background with 53% efficiency for the signal. A more detailed description of 
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the background suppression technique can be found in Ref. [9| and references therein. 

The understanding of the background that originates from other B meson decays is of 
great importance in the study of charmless B decays. We study the BB related back- 
ground using a large sample of MC generated BB generic events. We find that the 
dominant BB related background is due to B° — > D~tt + , D~ — > K^n" and due to 
B° -> J/^{2S))K% J/^{2S)) -> n+fj,- decays. We veto B° -> L>-7r+ events by 
requiring |M(iT^7r) — M D | > 100 MeV/c 2 . Modes with J/ip(ip{2S)) contribute due to muon- 
pion misidentification; the contribution from the J /ip(ip(2S)) — > e + e _ submode is found to 
be negligible after the electron veto requirement. We exclude J/il)(ip(2S)) background by 
requiring \M(ir + Tt-) li + li --Mj/i i \ > 70 MeV/c 2 and \M(7r + Tr~) li + lJ ,- —M^S)\ > 50 MeV/c 2 , 
with a muon mass assignment used here for the pion candidates. To suppress the background 
due to K/tt misidentification, we also exclude candidates if the invariant mass of any pair 
of oppositely charged tracks from the B candidate is consistent with the D + — > K$K + 
hypothesis within 15 MeV/c 2 (~ 2.5a), regardless of the particle identification information. 
The most significant background from charmless B decays is found to originate from the 
decay B° —>■ rfKg followed by v( — ► 7r + 7r~7. There is also a contribution from the two-body 
charmless decay B — > KgTr . Although this background is shifted by about 0.2 GeV from 
the AE signal region, it is important to take it into account for correct estimation of the 
background from continuum events. 

THREE-BODY SIGNAL YIELDS 

The AE distributions for B° — > Kg7T + TT~ candidates that pass all the selection require- 
ments are shown in Fig. ^a), where clear peak in the signal region is observed. The two 
dimensional AE versus Mb c plot is shown in Fig. IHb). In the fit to the AE distribution 
we fix the shape and normalization of the charmless BB background components from their 
measured branching fractions |10J and known number of produced BB events. For the BB 
generic component we fix only the shape and let the normalization float. The slope and nor- 
malization of the qq background component are free parameters. For signal we fix the width 
of the second Gaussian function at 31.0 MeV and the fraction at 0.19 as determined from 
MC simulation. The width of the main Gaussian is floating. The fit finds 1229 ± 62 signal 
events. The sigma of the main Gaussian is 15.3 ± 0.9 MeV. Results of the fit are shown in 
Fig. I^a), where different components of the background are shown separately for compari- 
son. There is a large increase in the level of the BB related background in AE < —0.15 GeV 
region. This is mainly due to B — > Dtt, D — > Kim decay. This decay mode produces the 
same final state as the studied process plus one extra pion that is not included in the energy 
difference calculation. The semileptonic decays B — > D^n, D — > Kfiu^ also contribute due 
to muon-pion misidentification. The shape of the BB background is described well by MC 
simulation. 

To examine possible quasi-two-body intermediate states in the observed B° — > Kg7i + n^ 
signal, we analyze the two-particle invariant mass spectra. To do so we define the B signal 
and sideband regions as illustrated in Fig. HJc). Defined in this way, the M^ c — AE side- 
bands are equivalent to the following sidebands in terms of the three-particle invariant mass 
M(Kit7t) and three-particle momentum P(Ktt7t): 

0.05 GeV/c 2 < \M(Ktttt) - M B \ < 0.10 GeV/c 2 ; P(Kirit) < 0.48 GeV/c 
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FIG. 1: AE distribution for the B° -» ^vr+vr" candidate events with \M hc - M B \ < 7.5 MeV/c 2 . 
Points with error bars are data; the upper curve is the fit result; the hatched histogram is the 
background, (b) Distribution of AE versus Mb c for the B° — > Kgir + Tr~ candidates in data, 
(c) Definition of the B signal and sideband regions in the M^ c — AE plane. 



and 

\M{Ktjtt) - M B \ < 0.10 GeV/c 2 ; 0.48 GeV/c < P(Ktttv) < 0.65 GeV/c. 

The B signal region is denned as an ellipse around the Mb c and AE mean values: 

(M hc -M B ) 2 AE 2 
(7.5 MeV/c 2 ) 2 (40 MeV) 2 

The efficiency of the requirements that define the signal region is 0.923. The total number 
of events in the signal region is 2207. The relative fraction of signal events in the B signal 
region is then determined to be 0.521 ± 0.025. 

The and tt + 7t~ invariant mass spectra for B° — * KgTT + 7r~ candidate events in the 

B signal region are shown as open histograms in Fig. |21 The hatched histograms show 
the corresponding spectra for background events in the M^ c — AE sidebands, normalized 
to the estimated number of background events. To suppress the feed-across between the 
7r + 7r~ and Kg7T^~ resonance states, we require the (tt + tt~) invariant mass to be larger 

than 1.5 GeV/c 2 when making the tc + tt~ (K^tt^) projection. The Kgii^ invariant mass 
spectrum is characterized by a narrow peak around 0.9 GeV/c 2 which is identified as the 
K* (892) ± and a broad enhancement around 1.4 GeV/c 2 . Possible candidates to assign to 
this enhancement are the scalar i^Q(1430) ± and tensor K^lAdQ)^ resonances. In the 7r + 7r~ 
invariant mass spectrum three distinct structures in the low mass region are observed. The 
most prominent one is slightly below 1.0 GeV/c 2 and is consistent with the /o(980). There 
is a clear indication for the p(770)° signal to the left of the /o(980) peak. Finally, there is 
a less prominent structure between 1.2 GeV/c 2 and 1.5 GeV/c 2 . We cannot identify unam- 
biguously the resonance state that is responsible for such a structure; possible candidates for 
a resonance state in this mass region mig ht be / (1370), / 2 (1270) and perhaps p(1450) H3- 
In what follows, we refer to this structure as /x(1300). It is worth noting that both Kir and 
mi two-body spectra in three-body B° — > Kg7r + 7r^ decays are similar to those observed in 
charged B meson decay to the three-body K + 7i + ti~ final state |1]. As we observe a clear 
B + — > XcoK + signal in the analysis of B + — > K + 7i + n~ decay, we expect B° — ► XcoK decay 
to occur at a similar rate. However, the 7r + n~ mass spectrum for the Xco region shown in 
Fig. 0(c) does not reveal the Xco signal clearly. 
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M(K°sjr) (GeV/c 2 ) M(7I + It ) (GeV/c 2 ) M(7I + It ) (GeV/c 2 ) 



FIG. 2: Two-particle invariant mass spectra for the B — ► i^^7r + 7r _ candidates in the -B signal 
region (open histograms) and for background events in the AE — Mt, c sidebands (hatched his- 
tograms), (a) M{K^ r K^ 1 ) spectrum with M(n + 7r~) > 1.5 GeV/c 2 (note that there are two entries 
per candidate in this plot); (b) M(vr + 7r") with M{K^ir c ) > 1.5 GeV/c 2 and (c) M(vr + 7r") in the 
XcO mass region with M(Koir ) > 1.5 GeV/c 2 . 



From these qualitative considerations it is apparent that an amplitude analysis is required 
for a more complete understanding of the individual quasi-two-body channels that contribute 
to the observed three-body B° — > Kg7i + ir~ signal. 



AMPLITUDE ANALYSIS 

In the preceding section we found that a significant fraction of the B° — > Kg7i + 7i^ signal 
can be assigned to quasi-two-body intermediate states. These resonances will cause a non- 
uniform distribution of events in phase space that can be analyzed using the technique 



pioneered by Dalitz . Multiple resonances that occur nearby in phase space will interfere 
providing an opportunity to measure their amplitudes and relative phases. This in turn 
allows us to deduce their relative fractions. 

The amplitude analysis of B meson three-body decays reported here is performed by 
means of an unbinned maximum likelihood fit. As the unbinned maximum likelihood fitting 
method does not provide a direct way to estimate the quality of the fit, we need a measure 
to assess how well any given fit represents the data. To do so the following procedure is 
applied. We first subdivide the entire Dalitz plot into 1 (GeV/c 2 ) 2 xl (GeV/c 2 ) 2 bins. If 
the number of events in the bin is smaller than N min = 16 it is combined with the adjacent 
bins until the number of events exceeds N min . After completing this procedure, the entire 
Dalitz plot is divided into a set of bins of varying size, and a x 2 variable for the multinomial 
distribution can be calculated as 

X 2 = -2 N iZnM{-)^ (1) 



=i 



71; 



where is the number of events observed in the i-th bin, and pi is the number of predicted 
events from the fit. For a large number of events this formulation becomes equivalent to 
the usual one. Since we are minimizing the unbinned likelihood function, our "x 2 " variable 
does not asymptotically follow a \ 2 distribution but it is bounded by a \ 2 variable with 
(Nuns — 1) degrees of freedom and a x 2 variable with (iVbins —k — 1) degrees of freedom 12l 



where k is the number of fit parameters. Because it is bounded by two x 2 variables, it should 
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FIG. 3: Dalitz plots for K^-k+tt candidates in the (a) AE — Mj-, c sidebands and (b) B signal 
region. 

be a useful statistic for comparing the relative goodness of fits for different models. A more 
detailed description of the technique is given in Ref. 



Fitting the Background Shape 

Before fitting the Dalitz plot for events in the signal region, we need to determine the 
distribution of background events. The background density function is determined from 
an unbinned likelihood fit to the events in the M^ c — AE sidebands defined in Fig. H^c). 
Figure Ol^a) shows Dalitz distribution for 8159 sideband events. This is about seven times 
the estimated number of background events in the B signal region. 

We use the following empirical parametrization to describe the distribution of background 
events over the Dalitz plot in the Kgir + TT~ final state 

B(s 13 ,s 23 ) = ai (e-P 1S12 +e" /3lS13 ) + a 2 e-^ S23 

+ 7l (\BW (K* (892)-)\ 2 + \BW(K* (892) + )\ 2 ) + 72 | J Biy(p(770)°)| 2 , (2) 

where s 12 = M 2 (K%ir~), s 13 = M 2 (K%ir + ), s 23 = M 2 (n+n-) and ot< (a x = 1.0), A and ^ 
are fit parameters; BW is a Breit-Wigner function. The first three terms in Eq. (0) are 
introduced to describe the background enhancement in the two-particle low invariant mass 
regions. This enhancement originates mainly from e + e~ — > qq continuum events. Due to the 
jet-like structure of this background, all three particles in a three-body combination have 
almost collinear momenta. Hence, the invariant mass of at least one pair of particles is in 
the low mass region. In addition, it is often the case that two high momentum particles 
are combined with a low momentum particle to form a B candidate. In this case there are 
two pairs with low invariant masses and one pair with high invariant mass resulting in even 
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FIG. 4: Results of the best fit to the Kgir + ir~ events in the AE — Mb c sidebands shown as 
projections onto two-particle invariant mass variables. Points with error bars are data; histograms 
are fit results. 

stronger enhancement of the background in the corners of the Dalitz plot. This is taken into 
account by terms 4 — 6 in Eq. To account for the contribution from real K*(892) ± and 
p(770)° mesons, we introduce two more terms in Eq. (J2j), that are (non-interfering ) sq uared 
Breit-Wigner amplitudes, with masses and widths fixed at world average values (lOj. The 
projections of the data and fits for the background events are shown in Figs.E] The xV^bins 
value of the fit is 241.7/197. 

Fitting the Signal 

The Dalitz plot for K^tt + h~ events in the signal region is shown in Fig. E[b). There are 
2207 events in the signal region that satisfy all the selection requirements. In an attempt 
to describe all the features of the K^Tt^ and 7t + tt~ mass spectra mentioned above, we use 
matrix element (refered to as model Ktttt—Cj) similar to those constructed in the analysis 
of B + -> K+Ti+Ti- decay Q: 

M Cj (K° s n + n-) = a K *e iSK *A 1 (n + K° s 7r-\K*(892) + ) + q k J Sk o A (ir + K° s ir-\K*(U30) + ) 
+ a p e^M 1 (K°7r + 7r-|p(770) ) + a fo e iS foA Flatte (K s 7r + TT-\f (980)) 
+ a fx e lS fxAj(K° s n + 7r~\f x (1300)) + a Xc0 e lS ™ A (K° s n+ \ Xc0 ) 
+ A ni (K s ir + 7i-), (3) 

where the subscript J denotes the unknown spin of the /x(1300) resonance; amplitudes 
di, relative phases 5i, are fit parameters. The masses and widths of all resonances are 
fixed at either their world average values 0] or at values determined from the analysis 
of B + -> K+tt+ix- decay (/ (980) and / x (1300)). The /o(98p_) is parametrized with a 
coupled channel Breit-Wigner function (Flatte parametrization [13]). For the non-resonant 
amplitude A nT we use an empirical parametrization 

A nr (K° s ir + n-) = < r e- QSl3 e^ r + a° r e - QS23 e^\ (4) 

Finally, as we currently do not perform the flavor analysis of the other B meson, we cannot 
distinguish whether a B or B meson decays to the K°7r + 7r~ final state. Thus the signal 
PDF is a non- coherent sum 

ScAKh+v-) = \M Cj (K°s^-)\ 2 + |^c J (^7T- 7 r + )| 2 . (5) 
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FIG. 5: Results of the fit to Kgir + Tr~ events in the signal region with the model K tttt— Co- Points 
with error bars are data, the open histograms are the fit result and hatched histograms are the 
background components. Insets in (a) and (b) show the K*(892) — ^(1430) mass region in 
20 MeV/c 2 bins; inset in (c) shows the XcO mass region in 25 MeV/c 2 bins. 




cos(0 H ) cos(0 H ) 




FIG. 6: Helicity angle distributions for Kgir + Tr events in different regions: 

(a) K*(892)+ (0.82 GeV/c 2 < M^tt*) < 0.97 GeV/c 2 ); 

(b) K *(1430)+ (1.0 GeV/c 2 < M(K+ir-) < 1.76 GeV/c 2 ); 

(c) p(770)° (M(7r+7r-) < 0.90 GeV/c 2 ); 

(d) / (980) (0.90 GeV/c 2 < M(tt + tt-) < 1.06 GeV/c 2 ). Points with error bars are data, the open 
histogram is the fit result and the hatched histogram is the background component. Note that in 
plots (a) and (b) there are two entries per B candidate. 
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TABLE I: Results of the fit to K^ir + ir~ events in the B signal region with model Ktttt— Cq. The 
first quoted error is statistical and the second is the model dependent uncertainty. 







Parameter 




Mode 


Fraction, % 


Phase, ° 


Mass, GeV/c 2 


Width, GeV/c 2 


K*(892)+7T- 


11.8 ± 1.4±±g;| 


(fi-x-prH 

VJ I 11ACU J 


891 66 ('fi-ic-pd'l h 01 


0.0508 (fixed) [10] 


^(1430)+7r~ 


64.8 ± 3.9lg:3 


45 ± 9 + ?, 

— lO 


1.412 (fixed) [1Q] 


0.294 (fixed) [1£] 


^(770)°^° 


12.9 ± l.9tH 


-9 ± 28±f| 


0.7758 (fixed) [10] 


0.1503 (fixed) [10] 


/ (980)^° 


16.0 ± 3.4±^;2 


36 ± 34±|? 


0.950 (fixed) [2] 


5 ^ = 0.23 (fixed) [2] 
Sjcx = 0.73 (fixed) [2J 


XcoK 


n .0+0.42+0.02 
u - 4<5 -0. 17-0.06 


-99 ± 37±| 


3.415 (fixed) [1Q] 


0.011 (fixed) fl£] 




3.68 ± 2.16ig:§ 


-135 ±25tjt 


1.449 (fixed) [2] 


0.126 (fixed) [2] 


Non-Resonant 


41.9*5.1^5 


5f = -22 ± 8± 6 6 
5f = 175±30t^ 






Charmless Only 


99.3 ±0.4 ±0.1 









While fitting the data, we choose the K*(892) + tt~ signal as our reference by fixing its 
amplitude and phase (ccr* = 1 and 8k* = 0). The numerical values of the fit parameters 
are given in Table |TJ Figure El shows the fit projections and the data. In addition, helicity 
angle distributions for several regions are shown in Fig. El The helicity angle is defined as 
the angle between the direction of flight of the tt + (tt~) in the Kgir + (Kgir~) rest frame and 
the direction of Kgir system in the B rest frame. For the 7r + 7r~ system the helicity angle is 
denned with respect to the positively charged pion. Gaps visible in Figs. Eland El are due to 
vetoes applied on invariant masses of two-particle combinations. All plots shown in Figs. El 
and El demonstrate good agreement between data and the fit. We also try to fit the data 
assuming /x(1300) is a vector (tensor) state. In this case we ascribe mass and width of 
p(1450) (^2(1270)) from PDG [10] to it. If parametrized by a single resonant state the best 
fit is achieved with a scalar assumption. 

MODEL & SYSTEMATIC UNCERTAINTIES 

To estimate the model dependent uncertainty and test for the contribution of other possi- 
ble quasi-two-body intermediate states such as fT*(1410) + 7r~, if *(1680) + 7r~, K\ (1430) + tt~ 
or f 2 (1270)K° , we include an additional amplitude of either of these channels to model 
Ktttt— Co and repeat the fit to data. For none of these channels is a statistically significant 
signal found. We also use several alternative (though also empirical) parametrizations of 
the non-resonant amplitude to estimate the related uncertainty 

• A nr (K° s TT + TT-) = aYe- asis e iS " r ; 

• A m (K° s TT + TT-) = aYe- asi3 e i5 " r + a™e- aS2:i e iS " r + afe- asi2 e i5 s r ; 
. A m {K° s TT+TT-) = §- e^ r ±f 

The dominant sources of systematic error are listed in Table ITTl For the branching fraction 
of the three-body B° — > KgTT + TT~ decay, we estimate the systematic uncertainty due to 
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TABLE II: List of systematic errors (in percent) for the three-body B° — > K^ir^n branching 
fraction. 



OUUI C(J 


Hjl I Ul 




9 


PID 


2.0 


Kg reconstruction 


3.0 


Event Shape requirements 


2.5 


Signal yield extraction 


5.4 


Model 


2.2 


MC statistics 


1.0 


Luminosity measurement 


1.0 


Total 


7.7 



possible losses from the mass cuts used to remove contributions from charmed particles by 
varying the relative phases and amplitudes of the quasi-two-body states within their errors. 
The systematic uncertainty due to requirements on event shape variables is estimated from 
a comparison of the | cos# t hr| and T distributions for signal MC events and B + — ► D°7r + 
events in the data. The uncertainty from the particle identification efficiency is estimated 
using pure samples of kaons and pions from the D° — > K~ir + decays, where the D° flavor is 
tagged using D* + — ► D°n + decays. We estimate the uncertainty due to the signal AE shape 
parameterization by varying the parameters of the fitting function within their errors. The 
uncertainty in the background parameterization is estimated by varying the relative fraction 
of the BB background component and the slope of the qq background function within their 
errors. The overall systematic uncertainty for the three-body branching fraction is estimated 
to be ±7.7%. 



RESULTS 

In previous sections we determined the relative fractions of various quasi-two-body in- 
termediate states in the three-body B° — > Kg7r + ir~ decay. To translate these numbers into 
absolute branching fractions, we first need to determine the branching fractions for the three 
body decay. To determine the reconstruction efficiency for the B° — > Kg7i + 7i^ decay, we use 
MC simulation where events are distributed over the phase space according to the matrix 
elements of the model Knir— Co- The corresponding reconstruction efficiency is 6.71±0.03% 
(including K° — > 7r + 7r~ fraction). 

Results of the branching fraction calculations are summarized in Table ILnl For final states 
where no statistically significant signal is observed we calculate 90% confidence level upper 
limits /go for their fractions. To determine the upper limit we use the following formula 

f/ 90 G(a, s: x)dx , , 

0-90 = J ° ) ' . , 6 

where G(a, s; x) is a Gaussian function with mean a and sigma s which are the measured 
mean value for the signal fraction and its statistical error. To account for the systematic 
uncertainty we decrease the reconstruction efficiency by one standard deviation. 
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TABLE III: Summary of branching fraction results. The first quoted error is statistical, the second 
is systematic and the third is the model error. 



Mode B(B -> Rh) x B(R -» hh) x 10 6 B{B -> ita) x 10' 



iig7r~*~7T charmless total 






4f .0 dz z.4 ±0.1 


A (oyzj^7T , a (^oyzj^ — » 


ii 7T 


O.Oi ± U. 11 ± U.4o_q 


o.4z zb l.Uo ± U.OO_q 43 


iig(1430)^7r , Aq(1430)^ - 


7^0 4- 
-> K TC 


O 1 C~\ A 1 O /I -r~U.o 

30.8 ± 2.4 ± 2.4Z3 


49.7 ± 3.8 ± 3.8_ 4 ; 8 


A"*(1410) + 7r-, K*(UW)+ - 




< 3.8 




K*(1680)+ir-, fl"*(1680)+ - 




< 2.6 




i^ 2 *(1430)+^-, i^ 2 *(1430)+ - 




< 2.1 




p(770)°K°, /5(770)° -» vr+vr 




6.13 ± 0.95 ± 0.47±i:8g 


6.13 ± 0.95 ± 0.47±i]gg 


/ (980)^°, / (980) -> tt+tt 




7.60 ± 1.66 ± 0.59±g;|f 




/ 2 (1270)K°, / 2 (1270) -> Tr"* 




< 1.4 




Non-resonant 






19.9 ± 2.5 ± 1.5i?;£ 






< 0.56 


< 113 



DISCUSSION & CONCLUSION 

With a 357 fb _1 data sample collected with the Belle detector, an amplitude (Dalitz) 
analysis of B meson decays to three-body charmless Kg7T + 7T~ final states is performed for 
the first time. Clear signals are observed in the B° — > A*(892) + 7r _ , B° — > Kq (1430) + 7r~, 
B° — > p(770)°A' and B° — > f (980)K° decay channels. The model uncertainty for these 
channels is small due to their narrow width and, for vector-pseudoscalar decays, due to clear 
signature because of the 100% polarization of the vector meson. Among these quasi-two- 
body channels the decay B° — ► p(770)°K° is observed for the first time. 

The branching fraction measured for the decay B° — > K*(892) + tt~ is in agreement with 
results from the amplitude analysis of the three-body B° — ► K + tt~'k decay [3,0], where 
the K*(892) + is reconstructed in the K + tt° channel. 

We also check possible contributions from B° — > l£j:(1430) + 7r-, B° — > fr*(1410) + 7r~, 
B° -> A"*(1680) + 7i- and B° -> / 2 (1270)i^° decays. We find no statistically significant 
signal in any of these channels and set 90% confidence level upper limits for their branching 
fractions. In the factorization approximation, charmless B decays to final states with a 
tensor state are expected to occur at the level of ~ 10~ 7 

We cannot identify unambiguously the broad structures observed in the M(7r + 7r~) ~ 
1.3 GeV/c 2 mass region. If approximated by a single resonant state, the best description 
is achieved with a scalar amplitude whose mass and width are consistent with those for 
the /o(1370). Because of the large uncertainty in / (1370) parameters and its coupling to 
7r + 7r~ [10|, a more definite conclusion is not possible at present. 
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